Force-induced unzipping of DNA in the presence of solvent molecules.

The melting of double-stranded DNA (dsDNA) in the presence of solvent molecules is a fundamental process with significant implications for understanding the thermal and mechanical behavior of DNA and its interactions with the surrounding environment. The solvents play an essential role in the structural transformation of DNA subjected to a pulling force. In this study, we simulate the thermal and force induced denaturation of dsDNA and elucidate the solvent dependent melting behavior, identifying key factors that influence the stability of DNA melting in presence of solvent molecules. Using a statistical model, we first find the melting profile of short heterogeneous DNA molecules in the presence of solvent molecules in Force ensemble. We also investigate the effect of solvent's strengths on the melting profile of DNA. In the force ensemble, we consider two homogeneous DNA chains and apply the force on different locations along the chain in the presence of solvent molecules. Different pathways manifest the melting of the molecule in both ensembles, and we found several interesting features of melting DNA in a constant force ensemble, such as lower critical force when the chain is pulled from the base pair close to a solvent molecule. The results provide new insights into the force-induced unzipping of DNA and could be used to develop new methods for controlling the unzipping process. By providing a better understanding of melting and unzipping of dsDNA in the presence of solvent molecules, this study provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA nanostructures.

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence.

Double-stranded RNA is the end-product of template-based replication, and is also the functional state of some biological RNAs. Similarly to proteins and DNA, they can be denatured by temperature, with important physiological and technological implications. Here, we use an in silico strategy to probe the thermal denaturation of RNA duplexes. Following previous results that were obtained on a few different duplexes, and which nuanced the canonical 2-state picture of nucleic acid denaturation, we here specifically address three different aspects that greatly improve our description of the temperature-induced dsRNA separation. First, we investigate the effect of the spatial distribution of weak and strong base-pairs among the duplex sequence. We show that the deviations from the two-state dehybridization mechanism are more pronounced when a strong core is flanked with weak extremities, while duplexes with a weak core but strong extremities exhibit a two-state behavior, which can be explained by the key role played by base fraying. This was later verified by generating artificial hairpin or circular states containing one or two locked duplex extremities, which results in an important reinforcement of the entire HB structure of the duplex and higher melting temperatures. Finally, we demonstrate that our results are little sensitive to the employed combination of RNA and water forcefields. The trends in thermal stability among the different sequences as well as the observed unfolding mechanisms (and the deviations from a two-state scenario) remain the same regardless of the employed atomistic models. However, our study points to possible limitations of recent reparametrizations of the Amber RNA forcefield, which sometimes results in duplexes that readily denature under ambient conditions, in contradiction with available experimental results.

Properties of phosphoramide benzoazole oligonucleotides (PABAOs). I. Structure and hybridization efficiency of N-benzimidazole derivatives.

In this work, we for the first time conducted a detailed study on the structure, dynamics, and hybridization properties of N-benzimidazole group-bearing phosphoramide benzoazole oligonucleotides (PABAOs) that we developed recently. By circular dichroism we established that the introduction of the modifications does not disrupt the B conformation of the DNA double helix. The formation of complexes is approximated by a two-state model. Complexes of PABAOs with native oligodeoxriboynucleotides form efficiently, and the introduction of such modifications reduces thermal stability of short duplexes (8-10 bp) by ∼5°Ð¡ per modification. Using UV-spectroscopy analysis, a neutral charge of the phosphate residue modified by the N-benzimidazole moiety in the pH range of 3-9.5 was found. The results confirm possible usefulness of PABAOs for both basic research and biomedical applications.

The binding of a c-MYC promoter G-quadruplex to neurotransmitters: An analysis of G-quadruplex stabilization using DNA melting, fluorescence spectroscopy, surface-enhanced Raman scattering and molecular docking.

The interactions of several neurotransmitter and neural hormone molecules with the c-MYC G-quadruplex DNA sequence were analyzed using a combination of spectroscopic and computational techniques. The interactions between indole, catecholamine, and amino acid neurotransmitters and DNA sequences could potentially add to the understanding of the role of G-quadruplex structures play in various diseases. Also, the interaction of the DNA sequence derived from the nuclear hypersensitivity element (NHE) III1 region of c-MYC oncogene (Pu22), 5'-TGAGGGTGGGTAGGGTGGGTAA-3', has added significance in that these molecules may promote or inhibit the formation of G-quadruplex DNA which could lead to the development of promising drugs for anticancer therapy. The results showed that these molecules did not disrupt G-quadruplex formation even in the absence of quadruplex-stabilizing cations. There was also evidence of concentration-dependent binding and high binding affinities based on the Stern-Volmer model, and thermodynamically favorable interactions in the form of hydrogen-bonding and interactions involving the π system of the aromatic neurotransmitters.

Equilibrium melting probabilities of a DNA molecule with a defect: An exact solution of the Poland-Scheraga model.

In this study we derive analytically the equilibrium melting probabilities for basepairs of a DNA molecule with a defect site. We assume that the defect is characterized by a change in the Watson-Crick basepair energy of the defect basepair, and in the associated two stacking energies for the defect, as compared to the remaining parts of the DNA. The defect site could, for instance, occur due to DNA basepair mismatching, cross-linking, or by the chemical modifications when attaching fluorescent labels, such as fluorescent-quencher pairs, to DNA. Our exact solution of the Poland-Scheraga model for DNA melting provides the probability that the labeled basepair, and its neighbors, are open at different temperatures. Our work is of direct importance, for instance, for studies where fluorophore-quencher pairs are used for studying single basepair fluctuations of designed DNA molecules.

Aberrant methylation scanning by quantitative DNA melting analysis with hybridization probes as exemplified by liquid biopsy of SEPT9 and HIST1H4F in colorectal cancer.

OBJECTIVE: The generally accepted method of quantifying hypermethylated DNA by qPCR using methylation-specific primers has the risk of underestimating DNA methylation and requires data normalization. This makes the analysis complicated and less reliable. METHODS: The end-point PCR method, called qDMA-HP (for quantitative DNA Melting Analysis with hybridization probes), which excludes the normalization procedure, is multiplexed and quantitative, has been proposed. qDMA-HP is characterized by the following features: (i) asymmetric PCR with methylation-independent primers; (ii) fluorescent dual-labeled, self-quenched probes (commonly known as TaqMan probes) covering several interrogated CpGs; (iii) post-PCR melting analysis of amplicon/probe hybrids; (iv) quantitation of unmethylated and methylated DNA alleles by measuring the areas under the corresponding melt peaks. RESULTS: qDMA-HP was tested in liquid biopsy of colorectal cancer by evaluating SEPT9 and HIST1H4F methylations simultaneously in the single-tube reaction. Differences in the methylation levels in healthy donors versus cancer patients were statistically significant (p < 0.0001), AUCROC values were 0.795-0.921 for various marker combinations. CONCLUSIONS: This proof-of-concept study shows that qDMA-HP is a simple, normalization-independent, quantitative, multiplex and "closed tube" method easily adapted to clinical settings. It is demonstrated, for the first time, that HIST1H4F is a perspective marker for liquid biopsy of colorectal cancer.

Thermodynamic Characterization of Nucleic Acid Nanoparticles Hybridization by UV Melting.

The advances in nucleic acid nanotechnology have given rise to various elegantly designed structural complexes fabricated from DNA, RNA, chemically modified RNA strands, and their mixtures. The structural properties of NA nanoparticles (NANP) generally dictate and significantly impact biological function; and thus, it is critical to extract information regarding relative stabilities of the different structural forms. The adequate stability assessment requires knowledge of thermodynamic parameters that can be empirically derived using conventional UV-melting technique. The focus of this chapter is to describe methodology to evaluate thermodynamic data of NANPs complexation based on DNA 12 base-pair (bp) duplex formation as an example.

Force induced DNA melting in the presence of an attractive surface.

The self avoiding walk (SAW) model of the polymer has been extended to study the equilibrium properties of double stranded DNA (dsDNA) where two strands of the dsDNA are modeled by two mutually attracting self-avoiding walks (MASAWs) in the presence of an attractive surface. We study simultaneous adsorption and force induced melting transitions and explore different phases of DNA. It is observed that melting is entropically dominated, which can be substantially reduced under the application of an applied force. We consider three scenarios, where the surface is weakly, moderately and highly attractive. For both weakly and moderately attractive surfaces, the DNA desorbs from the surface in a zipped form and acquires the conformation of a melted state with the rise in temperature. However, for a strongly attractive surface, the force applied at one end of the strand (strand-II) results in unzipping, while the other strand (strand-I) remains adsorbed on the surface. We identify this as adsorption-induced unzipping, where the force applied on a single strand (strand-II) can unzip the dsDNA if the surface interaction energy exceeds a specific threshold. We also note that at a moderate surface attraction, the desorbed-zipped DNA melts with an increase in temperature and the free strand (strand-I) gets re-adsorbed onto the surface.

Unlicensed origin DNA melting by MCV and SV40 polyomavirus LT proteins is independent of ATP-dependent helicase activity.

Cellular eukaryotic replication initiation helicases are first loaded as head-to-head double hexamers on double-stranded (ds) DNA origins and then initiate S-phase DNA melting during licensed (once per cell cycle) replication. Merkel cell polyomavirus (MCV) large T (LT) helicase oncoprotein similarly binds and melts its own 98-bp origin but replicates multiple times in a single cell cycle. To examine the actions of this unlicensed viral helicase, we quantitated multimerization of MCV LT molecules as they assembled on MCV DNA origins using real-time single-molecule microscopy. MCV LT formed highly stable double hexamers having 17-fold longer mean lifetime (τ, >1,500 s) on DNA than single hexamers. Unexpectedly, partial MCV LT assembly without double-hexamer formation was sufficient to melt origin dsDNA as measured by RAD51, RPA70, or S1 nuclease cobinding. DNA melting also occurred with truncated MCV LT proteins lacking the helicase domain, but was lost from a protein without the multimerization domain that could bind only as a monomer to DNA. SV40 polyomavirus LT also multimerized to the MCV origin without forming a functional hexamer but still melted origin DNA. MCV origin melting did not require ATP hydrolysis and occurred for both MCV and SV40 LT proteins using the nonhydrolyzable ATP analog, adenylyl-imidodiphosphate (AMP-PNP). LT double hexamers formed in AMP-PNP, and melted DNA, consistent with direct LT hexamer assembly around single-stranded (ss) DNA without the energy-dependent dsDNA-to-ssDNA melting and remodeling steps used by cellular helicases. These results indicate that LT multimerization rather than helicase activity is required for origin DNA melting during unlicensed virus replication.

Toward a Molecular Mechanism of Complementary RNA Duplexes Denaturation.

RNA duplexes are relatively rare but play very important biological roles. As an end-product of template-based RNA replication, they also have key implications for hypothetical primitive forms of life. Unless they are specifically separated by enzymes, these duplexes denature upon a temperature increase. However, mechanistic and kinetic aspects of RNA (and DNA) duplex thermal denaturation remain unclear at the microscopic level. We propose an in silico strategy that probes the thermal denaturation of RNA duplexes and allows for an extensive conformational space exploration along a wide temperature range with atomistic precision. We show that this approach first accounts for the strong sequence and length dependence of the duplexes melting temperature, reproducing the trends seen in the experiments and predicted by nearest-neighbor models. The simulations are then instrumental at providing a molecular picture of the temperature-induced strand separation. The textbook canonical "all-or-nothing" two-state model, very much inspired by the protein folding mechanism, can be nuanced. We demonstrate that a temperature increase leads to significantly distorted but stable structures with extensive base-fraying at the extremities, and that the fully formed duplexes typically do not form around melting. The duplex separation therefore appears as much more gradual than commonly thought.

C-Banding of Plant Chromosomes.

C-banding visualizes regions of chromosomes containing constitutive heterochromatin. It creates distinct patterns along the chromosome length and allows precise chromosome identification if C-bands are present in sufficient numbers. It is performed on chromosome spreads generated from fixed material, usually root tips or anthers. While there are numerous lab-specific modifications, all methods share the same steps: acidic hydrolysis, DNA denaturation in strong bases (usually saturated aqueous solution of barium hydroxide), washes in saline solution, and staining in Giemsa-type stain in a phosphate buffer. The method can be used for a wide range of cytogenetic tasks, from karyotyping, meiotic chromosome pairing analyses, to large-scale screening and selection of specific chromosome constructs.

In-silico prediction of RT-qPCR-high resolution melting for broad detection of emaraviruses.

Twenty-four species of RNA viruses contain members infecting economically important crops that are classified within the genus Emaravirus, family Fimoviridae. There are at least two other non-classified species that may be added. Some of these viruses are spreading rapidly and cause economically important diseases on several crops, raising a need for a sensitive diagnostic technique for taxonomic and quarantine purposes. High-resolution melting (HRM) has shown to be reliable for the detection, discrimination, and diagnosis of several diseases of plants, animals, and humans. This research aimed to explore the ability to predict HRM outputs coupled to reverse transcription-quantitative polymerase chain reaction (RT-qPCR). To approach this goal a pair of degenerate genus-specific primers were designed for endpoint RT-PCR and RT-qPCR-HRM and the species in the genus Emaravirus were selected to framework the development of the assays. Both nucleic acid amplification methods were able to detect in-vitro several members of seven Emaravirus species with sensitivity up to one fg of cDNA. Specific parameters for in-silico prediction of the melting temperatures of each expected emaravirus amplicon are compared to the data obtained in-vitro. A very distinct isolate of the High Plains wheat mosaic virus was also detected. The high-resolution DNA melting curves of the RT-PCR products predicted in-silico using uMeltSM allowed saving time while designing and developing the RT-qPCR-HRM assay since the approach avoided extensive searching for optimal HRM assay regions and rounds of HRM tests in-vitro for optimization. The resultant assay provides sensitive detection and reliable diagnosis for potentially any emaravirus, including new species or strains.

Pulling short DNA with mismatch base pairs.

Due to misincorporation during gene replication, the accuracy of the gene expression is often compromised. This results in a mismatch or defective pair in the DNA molecule (James et al. 2016). Here, we present our study of the stability of DNA with defects in the thermal and force ensembles. We consider DNA with a different number of defects from 2to16 and study how the denaturation process differs in both ensembles. Using a statistical model, we calculate the melting point of the DNA chain in both the ensemble. Our findings display different manifestations of DNA denaturation in thermal and force ensembles. While the DNA with defects denatures at a lower temperature than the intact DNA, the point from which the DNA is pulled is important in force ensemble.

Choline ester based ionic liquid: A multi-functional system to enhance nucleic acid stability, drug solubilization and cell penetration.

Ionic liquids (ILs) are emerging systems with applications in varying areas of biomedical research. This study aims at developing a biocompatible, dual function choline ester-based IL with chloride as anion ([Ch] IL) for stabilizing nucleic acids (DNA) and enhancing cellular uptake of drugs. The ability of IL to complex with DNA was characterized using electrophoresis, dye displacement and UV absorbance. The effect of pH on complex stability and protection of DNA from nuclease were also studied. Even though [Ch] IL had positive zeta potential and showed effective complex formation, at physiological pH the zeta potential of the complex decreased and became negative, thereby, destabilizing the complex. To address this, citric acid (CA) was added to [Ch] IL which facilitated strong complexation. Further, DNA could be retrieved from these complexes without compromising its purity and integrity. Additionally, [Ch] IL was found to improve the cellular uptake of doxorubicin by improving its solubility in water. Thus, we demonstrate that the [Ch] IL developed here can enhance nucleic acid stability, drug solubilization and cell penetration. Our results show that the developed [Ch] IL can be used for long term storage of nucleic acids as well as for enhancing permeation of drugs in vivo.

Interaction of an anticancer benzopyrane derivative with DNA: Biophysical, biochemical, and molecular modeling studies.

BACKGROUND: SIMR1281 is a potent anticancer lead candidate with multi- target activity against several proteins; however, its mechanism of action at the molecular level is not fully understood. Revealing the mechanism and the origin of multitarget activity is important for the rational identification and optimization of multitarget drugs. METHODS: We have used a variety of biophysical (circular dichroism, isothermal titration calorimetry, viscosity, and UV DNA melting), biochemical (topoisomerase I & II assays) and computational (molecular docking and MD simulations) methods to study the interaction of SIMR1281 with duplex DNA structures. RESULTS: The biophysical results revealed that SIMR1281 binds to dsDNA via an intercalation-binding mode with an average binding constant of 3.1 × 106 M-1. This binding mode was confirmed by the topoisomerases' inhibition assays and molecular modeling simulations, which showed the intercalation of the benzopyrane moiety between DNA base pairs, while the remaining moieties (thiazole and phenyl rings) sit in the minor groove and interact with the flanking base pairs adjacent to the intercalation site. CONCLUSIONS: The DNA binding characteristics of SIMR1281, which can disrupt/inhibit DNA function as confirmed by the topoisomerases' inhibition assays, indicate that the observed multi-target activity might originate from ligand intervention at nucleic acids level rather than due to direct interactions with multiple biological targets at the protein level. GENERAL SIGNIFICANCE: The findings of this study could be helpful to guide future optimization of benzopyrane-based ligands for therapeutic purposes.

Sequence-specific binding behavior of coralyne toward triplex DNA: An ultrafast time-resolved fluorescence spectroscopy study.

Triplex DNA structure has potential therapeutic application in inhibiting the expression of genes involved in cancer and other diseases. As a DNA-targeting antitumor and antibiotic drug, coralyne shows a remarkable binding propensity to triplex over canonical duplex and thus can modulate the stability of triplex structure, providing a prospective gene targeting strategy. Much less is known, however, about coralyne-binding interactions with triplex. By combining multiple steady-state spectroscopy with ultrafast fluorescence spectroscopy, we have investigated the binding behaviors of coralyne with typical triplexes. Upon binding with a G-containing triplex, the fluorescence of coralyne is markedly quenched owing to the photoinduced electron transfer (PET) of coralyne with the G base. Systematic studies show that the PET rates are sensitive to the binding configuration and local microenvironment, from which the coexisting binding modes of monomeric (full and partial) intercalation and aggregate stacking along the sugar-phosphate backbone are distinguished and their respective contributions are determined. It shows that coralyne has preferences for monomeric intercalation within CGG triplex and pure TAT triplex, whereas CGC+ triplex adopts mainly backbone binding of coralyne aggregates due to charge repulsion, revealing the sequence-specific binding selectivity. The triplex-DNA-induced aggregation of coralyne could be used as a probe for recognizing the water content in local DNA structures. The strong π-π stacking of intercalated coralyne monomer with base-triplets plays an important role in stabilizing the triplex structure. These results provide mechanistic insights for understanding the remarkable propensity of coralyne in selective binding to triplex DNA and shed light on the prospective applications of coralyne-triplex targeted anti-gene therapeutics.

High-Resolution Melting (HRM) Genotyping.

High-resolution melting (HRM) analysis is a simple, fast, and inexpensive real-time polymerase chain reaction (PCR)-based method used to identify genetic variation between populations and detect single-nucleotide polymorphisms (SNPs) in nucleic acid sequences. HRM is a powerful technique that detects the differences between SNP allele melting temperatures by using a fluorescent dye inserted into the duplex deoxyribonucleic acid (DNA) structure. Prior to performing HRM analysis, optimizing the primer design, PCR mixture, and software settings is essential to obtain accurate and reliable results. In this chapter, we describe a detailed SNP genotyping method that includes primer design and the analysis of the shapes and positions of the melt curve of the luminescence intensity of the fluorescent dye attached to amplified DNA using software of qPCR instruments. This protocol is applicable for genotyping germplasm, genetic mapping, and marker-assisted breeding in plants.

Modified High-Resolution Melting (HRM) Marker Systems Increasing Discriminability Between Homozygous Alleles.

Targeted single-nucleotide polymorphism (SNP) genotyping, especially for functional nucleotide polymorphism, is widely used for current breeding programs in crops. One of the cost- and time-effective approaches for genotyping is high-resolution melting (HRM) analysis for polymerase chain reaction (PCR) amplicons, including target SNP. The reliability of a genotype obtained from an HRM marker depends on the difference in Tm values between two amplicons. Increasing the reliability of HRM marker genotypes could be archived with the selection of the best nearest neighboring nucleotide substitution (NNNs) in primer sequences surrounding SNPs. This chapter provides an easy-way protocol to design primer sequences for NNNs-HRM markers with table and web service, as well as several tips to develop HRM markers that distinguish between homozygous alleles (e.g., between A/A and C/C).

Dinuclear complex-induced DNA melting.

Dinuclear copper complexes have been designed for molecular recognition in order to selectively bind to two neighboring phosphate moieties in the backbone of double strand DNA. Associated biophysical, biochemical and cytotoxic effects on DNA were investigated in previous works, where atomic force microscopy (AFM) in ambient conditions turned out to be a particular valuable asset, since the complexes influence the macromechanical properties and configurations of the strands. To investigate and scrutinize these effects in more depth from a structural point of view, cutting-edge preparation methods and scanning force microscopy under ultra-high vacuum (UHV) conditions were employed to yield submolecular resolution images. DNA strand mechanics and interactions could be resolved on the single base pair level, including the amplified formation of melting bubbles. Even the interaction of singular complex molecules could be observed. To better assess the results, the appearance of treated DNA is also compared to the behavior of untreated DNA in UHV on different substrates. Finally, we present data from a statistical simulation reasoning about the nanomechanics of strand dissociation. This sort of quantitative experimental insights paralleled by statistical simulations impressively shade light on the rationale for strand dissociations of this novel DNA interaction process, that is an important nanomechanistic key and novel approach for the development of new chemotherapeutic agents.

A predictive model for the thermomechanical melting transition of double stranded DNA.

By extending the classical Peyrard-Bishop model, we are able to obtain a fully analytical description for the mechanical response of DNA under stretching at variable values of temperature, number of base pairs and intrachains and interchains bonds stiffness. In order to compare elasticity and temperature effects, we first analyze the system in the zero temperature mechanical limit, important to describe several experimental effects including possible hysteresis. We then analyze temperature effects in the framework of equilibrium Statistical Mechanics. In particular, we obtain an analytical expression for the temperature-dependent melting force and unzipping assigned displacement in the thermodynamical limit, also depending on the relative stability of intra vs. inter molecular bonds. Such results coincide with the purely mechanical model in the limit of zero temperature and with the denaturation temperature that we obtain with the classical transfer integral method. Based on our analytical results, we obtain explicitly phase diagrams and cooperativity parameters, where also discreteness effect can be accounted for. The obtained results are successfully applied in reproducing the thermomechanical experimental melting of DNA and the response of DNA hairpins. Due to the generality of the model, exemplified in the proposed analysis of both overstretching and unzipping experiments, we argue that the proposed approach can be extended to other thermomechanically induced molecular melting phenomena. STATEMENT OF SIGNIFICANCE: We obtain a fully analytical description of the complex wiggly energy landscape of two stranded macromolecules under unzipping loading. Based on Equilibrium Statistical Mechanics, we describe the combined thermomechanical effects and the melting transition of double stranded molecules such as nucleic acids. This is proved by quantitatively predicting the experimental behavior of both melting of DNA and DNA hairpins opening. While analytical results have been previously attained under special conditions on the relative stiffness of the covalent vs. non-covalent bonds of the base pairs, our model is completely general in this respect, thus representing a tool in the perspective of the design at the molecular scale. We show that the obtained model can be fully inscribed in the theory of phase transitions giving a new interpretation of the thermomechanical behavior of double stranded molecules.